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f ^ ■ Abstract 

>— ' . Proposals for nonlinear extenstions of quantum mechanics are dis- 

cussed. Two different concepts of "mixed state" for any nonlinear version 
of quantum theory are introduced: (i) genuine mixture corresponds to op- 
erational "mixing" of different ensembles, and (ii) a mixture described by 

rS—t ' single density matrix without having a canonical operational possibility 

C*H, to pick out its specific convex decomposition is called here an elementary 

mixture. Time evolution of a class of nonlinear extensions of quantum 
mechanics is introduced. Evolution of an elementary mixture cannot be 
generally given by evolutions of components of its arbitrary convex de- 
compositions. The theory is formulated in a "geometric form": It can 
be considered as a version of Hamiltonian mechanics on infinite dimen- 
sional space of density matrices. A quantum interpretation of the theory 
is sketched. 



1 Introduction 

In popular and well written book |jj on conceptual foundations of quantum me- 
chanics (QM) there is a subsection [Chap. 9-4, p. 278] on nonlinear Schrodinger 
equation inserted into the section entitled Some impossible processes, and con- 
taining a proof of inconsistency of any nonlinear evolution in QM "if we retain 
the other postulates of quantum theory without any change" . These "other pos- 
tulates" include statistical interpretation and time evolution of density matrices 
corresponding to some of our genuine mixtures, as well as the definition of en- 
tropy by von Neumann which is appropriate just to our elementary mixtures 

*The talk presented at the conference: "New Insights in Quantum Mechanics", Goslar, 
Aug. 31st - Sept. 4th, 1998. 



(cf. Subsection |2.l| ). The nonlinear extensions of QM appearing usually in lit- 
erature 1^, [| [| are formulated for time evolution of vector states (i.e. wave 
functions); trials to extend such dynamics to evolutions of density matrices led 
to inconsistencies, cf. e.g. criticism [|J of ||. 

We argue that (at least some of) the inconsistencies mentioned above are 
consequences of improper interpretation of "mixtures" and of definition of their 
time-evolution. We shall introduce two concepts of "mixed states" in this con- 
tribution, as well as a class of nonlinear extensions of dynamics avoiding the 
mentioned inconsistencies. Our formulation of nonlinear quantum theory |6[ 
was inspired by symplectic reformulation of QM M, and was encouraged by 
publication of the Weinberg's proposal ||. 

2 Two concepts of "mixture" and nonlinearity 
in QM 

In standard elementary formulation of QM (linear, without superselection rules, 
describing systems with "finite number of degrees of freedom" ) the set of bounded 
"observables" is described by the set C{TL) S of all bounded selfadjoint operators 
in a complex Hilbert space 7i, and the set of "states" is described by the set of 
density matrices <S* := T(Tt)+i consisting of positive trace-class operators of 
unit trace. Any density matrix g g 5* is either of the form 

Q = P (p =P* = P;e£(H) s , <peH, (l) 

corresponding to a "pure state" , or it is a nontrivial convex combination of one 
dimensional projections P v . Each density matrix corresponding to a nontrivial 
mixture, i.e. g 2 ^ g, can be written in infinitely many mutually different ways 
in the form of its convex decompositions: 




Q = 2_, \Qi ( = 2^ K k v k ) , Qj, Vk e T(H) + i, Xj,K k e [0, 1) c R+. (2) 

The density matrices of the form (0) represent fil "mixed states" , resp. "mix- 
tures" . 

Each such a state is uniquelly determined by determination of expectations 
of "sufficiently many" observables A G C(TL) S , i.e. by numbers 

(A) e :=Tr(g-A), VAe£(H) s . (3) 

Such a set of state-determining expectations corresponding to a nontrivial mix- 
ture (H) can arise in QM in two different experimental situations corre 



sponding to two possible interpretations of density matrices, cf. Subsection 2.1 
These two possibilities are not empirically distinguishable in the standard frame- 
work of QM, but their distinction becomes crucial, as we shall show in this paper, 
in any nonlinear generalization of the theory. 



2.1 Two concepts of mixed states 

Let us distinguish here two kinds of preparation procedures for states described 
in QM by a density matrix g, and let us also introduce the corresponding two 
concepts of "mixed states" . 

(i) . A state of a given system in QM described by a density matrix g ^ g 2 is of- 
ten interpreted as representation of the statistical ensemble {g} of copies of the 
considered physical system consisting of subensembles numbered by an index 
j G J(= an index set), each of which is prepared by its own macroscopicaly dis- 
tinguishable preparation procedure preparing the system in a state Qj G S* , and 
occuring in {g} with probability Aj, g = J^. \jgj- Hence the term "mixture" 
corresponds here to intuitive notion of "mixing" of a set of copies of the system 
occuring in states gj with relative frequencies Xj. There are infinitely many of 
such mutually different (possibly quantummechanicaly mutually incompatible) 
preparation procedures leading to the same density matrix g which are in the 
framework of QM indistinguishable Q by measurements of the quantities (||) . 
Each of the "components" Qj such that Qj ^ g 2 could be decomposed further 
into "less mixed" components, Qj = J^k^jkQjk, etc. Let the components Qj of 
the decomposition corresponding to the described preparation procedure of the 
ensemble {g} be empirically indecomposable, i.e. they are elementary mixtures 
in the sense of the next point (ii) . Then the state corresponding to the ensemble 
{g} will be called a genuine mixture. It is characterized by a specific convex 
combination of elementary mixtures: {g} = {g}\ := {(Xj;gj) : j G </}. Other 
such decompositions of the same density matrix will correspond to different 
genuine mixtures (empirically distinguishable in some generalizations of QM, 
cf. later in this paper). 

(ii). Let a density matrix g G <S» associated with the considered system S be 
given. It is always possible (at least mathematically) to find another QM-system 
S' and a vector state $ G Ti <8> H' of the combined system S + S' such that its 
restriction to the subsystem S is g: 

Tr(P$-A®I nl )=Tr(g-A), \/Ae£{H) s , (4) 

where In 1 £ C(Ti') is the unit observable of the attached system S' (called 
sometimes an "ancilla"). 

Let us assume that the state g of S is obtained by such a restriction from an 
empirically prepared state P$ of a system S+S', the observables of which contain 
all elements of C(7i <E> H') s (hence, no superselection rules). The systems S and 
S' are dynamically independent. The state g of S (and the corresponding statis- 
tical ensemble of copies of S) prepared in this way will be called an elementary 
mixture. These states are a special case of genuine mixtures corresponding to 
the trivial decomposition of g into a single element. 

This second method of preparation of a state g do not specify any decom- 
position of g into "simpler" , or "purer" states. We need not specify in the 
following a way of preparation of an elementary mixture; in the forthcoming 



considerations, any elementary mixture g will be a "full fledged" description of 
a quantum state, in an equal degree for nontrivial mixtures g 2 ^ g, as well as 
for pure states. From the point of view of physical intuition, the elementary 
mixtures are distinguished by impossibility to specify their decomposition to 
subensembles determined by their physical preparation. 

Remark 1 We shall see how nonlinear quantum dynamics can distinguish dif- 
ferent genuine mixtures with the same "bary centre" g by measuring the quanti- 
ties in ([Jj only. Detailed specification of conditions under which we are dealing 
either with "genuine", or with "elementary" mixture will be actual only after 
observing some nonlinear quantum evolution: It is also a question for experi- 
mentalists. Here we stress just theoretical possibilities and their connections. 

These different possible interpretations of density matrices are usually ignored 
because of their experimental indistinguishability in (linear) QM.R They could 
be distinguished by considering also states of preparation apparatuses, i.e. by 
considering correlations of states of S with some "macroscopic parameters" (we 
shall not go into details of this point here, cf., however, Remark g). 

Remark 2 Te above physically intuitive determination of difference between the 
two concepts of "mixed states " can be made mathematically clear after accepting 
an extended set of "observables" for the considered system S. Let the C* -algebra 
of observables be the setC s {S*,C{TL)) of continuous (in some conveniently spec- 
ified topologies) functions a : v ^^ &(f) G £(H) on "conveniently" compactified 
5*, 5* 3 v, with bounded operator values. The centre of this C* -algebra contains 
classical quantities consisting of the commutative C* -algebra of scalar-valued 
continuous functions a £ C(5*). The elementary mixtures are then just the pure 
states of this commutative subalgebra . Intuitively, this extension of the set of 
observables corresponds to the extension of the physical system S by macroscopic 
parameters of (preparation procedures of) its (microscopic) states. 

2.2 Nonlinear transformations distinguish the two kinds 
of mixed states 

The above mentioned indistinguishability of different decompositions (g) of any 
given mixture described by a density matrix g is conserved by linear time evo- 
lutions (or any unitary transformations) 

cpt(g) := u t £u£ , u t = exp(-itH). (5) 

We have then for all t £ R: 

Vt(o) = Yl ^oVtiQj) = Yl K Wt{vk), <Pt=o{p) = Q- (6) 

j k 



1 The only situations where the present author noticed el discussion of these different inter- 
pretations were works on quantum measurement problem [H. The genuine mixtures are called 
there also Gemenge, and the interpretation: ignorance interpretation. 



Conversely, if a continuous group of invertible transformations ipt of 5* satis- 
fies (^|) for all decompositions (^) of all density matrices g, then it is described 
by some unitary operators u* = exp(— iti/),[|| [Theorem 3.2.8, and Example 
3.2.14]. In these cases, determination of time evolution for pure states is suffi- 
cient to determine evolution of all states uniquely. 

Let us assume that ipt is a nonlinear family of transformations of <S*, i.e. (|6|) 
is no more identically valid. Hence there is g and its decompositions (||) such 
that for a "nonlinear t" : 

(Me) ¥>) ^^jPtiQj) ^ ^Kfe^(^)- (7) 



In other words, the mappings ipt : 5* i—> <S* are not all affine now, and the 
evolution of components gj of a decomposition (|2|) of a density matrix g need 
not determine a "corresponding" evolution of g. Evolutions of elementary mix- 
tures have to be specified on the whole space of quantum states iS». Moreover, 
eq. (ffl) shows that nonlinear evolution can distinguish between different genuine 
mixtures corresponding to the same initial density matrix. 

3 Nonlinear Extensions of Quantum Dynamics 

Let us specify here a class of evolutions generalizing the linear ones (ph. We 
shall ignore here problems with unboundcdness of generators. After defining 
an evolution of all elementary mixturesn we shall stress its distinction from the 
corresponding evolutions for genuine mixtures. We shall mention briefly the 
statistical interpretation of the generalized quantum theory. More details on 
the theory can be found in 0. 

3.1 Dynamics on the quantum phase space 

Let us consider the space <S* of all density matrices as the quantum phase 
space. The dynamics on 5* is just Hamiltonian (better: Poisson) classical 
dynamics on the infinite-dimensional submanifold (with boundary) 5* of the 
linear space T S (TL) (= symmetric trace class) endowed with the trace norm and 
with the Poisson brackets 

{f,h}(g)=iTr(g[D e f,D e h]), f,h € C7°°(T S (W),R). (8) 

Here D e f is the differential of / 6 C° c (T s (Ti.), R) considered (in a canonical 
way ||, cf. Remark | as a bounded linear operator on H, and [•, •] is the 
commutator in C{TL). 

2 Similar dynamics for density matrices was formulated by Czachor et al. [hd] in terms of 
Lie-Nambu brackets. 



Remark 3 The correspondence between D g f and operators in C(TC) S is given by 
the duality between the space of selfadjoint trace class operators T s (7i) (9 g, v), 
and its dual C{TC) (3 A, B) representing (real-)linear junctionals v t— > (B; v) G 
R on T S (TL) so that (B;v) :— Tr(B-v). Hence the differential D e f calculated 
on the vector v according to its definition as a linear functional on T S (H), can 
be represented by the operator D g f G C(7i) s determined by the relation 



D e f(u) = (D e f iV ):=± 



f{ Q + tv)=Tr{vb e f). (9) 



t=o 



Being this valid for all v G %(H.) 7 the operator D e f is determined unambigu- 
ously by (pi]. We shall use D = D in the following. 

Let Q G C°° (7, (7i) , R) be a Hamiltonian for our dynamics (for "more realistic" 
cases this function Q is not everywhere defined and it is unbounded) H. Then 
the Hamiltonian flow (pjf can be described || by unitary cocycle uq consisting 
of the set of unitary operators {uq (t,g) :ieR,ge5,} on H satisfying a version 
of the nonlinear Schrodinger equation: 

*37 U Q(*' B) =D g ( t) Q-u Q (t, g), u Q (Q,g) = 0; (10) 

and fulfilling the cocycle identity: 

u Q (t + s,g) =u Q (s,(pf(g))u Q (t, g), Vs,t,eR. (11) 

Here enters the nonlinearity via 

g(t) = <p?(g) := u Q (t, g)gu Q (t, g)* , g{0) := g. (12) 

The equation ( |l0| ) leads to the evolution equation written directly for g(t): 

ij t g{t) = [D g(t) Q,g{t)]. (13) 

Let us rewrite eq. ( |l0|) into a common form of nonlinear Schrodinger equation 
for wave functions ip(t) :— ug(t, Pip)tp G H (we set D^, = Dp^): 

ij t ^(t) = D m Q-i,(t). (14) 

It is seen from (|l2|), and from unitarity of uq(£, g) that the evolution tp® leaves 
spectral characteristics of all density matrices invariant. Hence it conserves also 
purity of states. If the evolution were known for vector states only, evolution 
for density matrices would remain undefined. In specific cases, it is possible [nG| 
to find a "natural extension" of the Hamiltonian Q from the subset P{Ti) C 5* 
of vector states to other parts of 5*; mathematical guides for such an exten- 
sion might be, e.g. continuity, symmetry, or some aesthetical, resp. formal 
considerations. 



Example 1 Let us take, e.g. Tt=L 2 (R n ) with (ip\(p) :— f ip(x)(f(x)d n X. Let us 
write density matrices g "in the x-representation" with a help of their operator 
kernels q(x, y): 

\(nl>\{x)= ( ' Q (x,y)^(y)d n y, iP^H. (15) 



Projection operators P^ have the kernels P^{x,y) = ||^|| ip(x)ip(y). Let the 
Hamiltonian function Q will be taken as the (unbounded) functional 

Q(P^) ■-Tr(P^-H ) + -^ f P^x, x) a+1 d n x, (16) 

with Ho some selfadjoint (linear) operator on L 2 (R n ), and a > O.n Let t i— > 
P^n\,ip(0) := ip be any differentable curve through P^ G P{7i), and let P^, G 
Tp^P(H) be its tangent vector expressed by an operator according to Remarks. 
Then the (unbounded, nonlinear) Hamiltonian D^Q can be expressed by: 



Tr(D^Q-P 4 ,) := - 



Q(Pm)> ( 17 ) 



what leads to the corresponding form of "nonlinear Schrodinger wave-equation" 
foript '■= ^(t) : 



dt vt 



(x) = [HoA}(x)+£\i>t(x)\ 2a Mx), 11^11 = 1- (18) 



One possible extension of this nonlinear dynamics to the whole space <S* is ob- 
tained by "the substitution g y-^ P^p ", i.e. by the choice of the Hamiltonian 

Q(g) := Tr(g-Ho) + -^— f g{x, x) a+1 d n x, (19) 

a + 1 J 

and the corresponding dynamics is then described by (pqj with 

D e Q(v) = Tr{v-H ) + e f g(x, x) a v(x, x) d n x. (20) 



We shall compare in the next subsection the evolutions of the mixed states 
described by the same initial density matrix (0) for the two distinguished inter- 
pretations. 

3.2 Genuine mixtures, their dynamics and interpretation 

Since genuine mixtures {g} :— {g}\ :— {(\j]gj) : j E J} are considered as 
probability distributions over the space of "elementary events" <S* (in the sense of 

3 We shall proceed, in the presented "Examples", in a heuristic way, by "plausible" formal 
manipulations; the necessary mathematical comments are omitted here. 



the classical Kolmogorov probability theory) with probabilities Xj concentrated 
in the points Qj € 5* , adequate description of genuine mixtures are probability 
measures M{g} on <£*■ Let us denote S e the Dirac probability measures: S g (A) = 
O () ^ A C 5,. The measure /i{ e j describing the state-decomposition (0) 
describes the integral of functions / e L 1 (S*,/j,{ g }): 

/"{,}(/) = I>M/) = I>M). (21) 

In this language, the elementary mixtures might be considered as the special 
case of genuine ones, in which the decomposition {g} of g according to (H) is 
trivial: Ai := 1, g\ := g, Xj = (Vj ^ 1). 

If there is given a Hamiltonian time evolution ip^ on >S* , the evolution of the 
state described by ^{ g \ can be described by the corresponding evolution of the 
distinquished decomposition: 

M := {{X j]Qj ) :jEJ}^ {g} t := {(A,; tpf( Qj )) : j G J}. (22) 

This evolution is described equivalently by the evolution of measures [it on 5* : 

(t; m) !-»■ Mt = M ° V-t- (23) 

Let us illustrate, by explicit calculation, the difference between time evolutions 
of the same initial density matrix considered in its two different interpretations. 

Example 2 Let us take the system with its "extended" dynamics from Ex- 
ample uL and let us fix a nontrivial mixture g of several vector states P^, . : 
q = ^jXjP^p-. Let us calculate the difference between the derivatives with re- 
spect to the time in t — of the two evolutions: (i) of the genuine mixture evo- 
lution ^2- \j<pj? (Pi/ij), and (ii) of the elementary mixture evolution ipj~(g). We 
shall calculate the right side of (UjJ for the two cases and take their difference. 
Let us write the kernel "in x -representation" of g as the convex combination of 
the vector-state kernels: 

g(x,y) =£>||^||-Vi(z)^(y)- (24) 

j 

The (symbolic) "kernel" of the Hamiltonian D e Q can be written: 

D g Q(x, y) = H (x, y) + eS(x - y)g(x, x) a . 

Here, S(-) is the Dirac distribution on R™. We have to express the difference 
A} e (x,y) between the kernels (in x-representation) of the operators 

~Y^Xi[D^ t )Q,P^\, and [D s{t) Q, g(t)}, 



what expresses the difference between time derivatives of "the same density ma- 
trix" g — X^Aj-^V in the two interpretations. The linear operator Hq does not 
contribute into this difference. The kernels of commutators entering into the 
calculation are (for all v S <S* ) of the form 

[D„Q,v](x,y) = [H ,u](x,y) + ev(x,y)(v(x 1 x) a -v(y,y) a ). 

We can (and we shall) take all ||^>j|| = 1. Let us denote 



X^Hx) :=\^^)\ 2a ~ [J2 X ^^\" 



Ao\(„\ ._ |„/,./V\|2o! _ 

Then the wanted difference at t = is 

& {e} {x,y):=&l e \x 1 y)=sY,^U^ J (y)(x\ e} {x)-x\ e} {y))- (25) 

3 

By proving that the operator A/ e \ is not identical zero for all {g}, we can prove 
nontrivial difference of the two time evolutions explicitly. This can be easily 
proved for Xi := 1 — A2, and ipi,ip2 specific two-valued functions concentrated 
on disjoint subsets o/R™. 

Remark 4 Our proposal for abstract interpretation scheme of the theory 
looks as follows: 

The set of observables is C S (S*, £(Tt)) (cf Remark^). Their expectations in 
the states /i{ e j are 

(&> {c} s ^AjTrG&aUfc)). (26) 

3 

The higher momenta are (a fe ){ e }. k = 2, 3, . . ., what allows us to calculate prob- 
ability distributions of (also microscopic) observables. The "classical" distribu- 
tion of elementary mixtures composing a genuine one can be given by an arbi- 
trary (not only discrete) measure \i on 5* (endowed by a Borel structure W). 
Then we have 

(&>„= / Tr( ga (g))fi(dg). (27) 

Js, 

The system is here described as a composition of a QM-system on TL, and a 

classical system with phase space <S*. The above considered nonlinear dynamics 

of the quantum system can be described as linear dynamics of this composed 

system |g/. 
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